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Abstract 

We analyze theoretically and experimentally the electronic structure and charging diagram of 
three coupled lateral quantum dots filled with electrons. Using the Hubbard model and real-space 
exact diagonalization techniques we show that the electronic properties of this artificial molecule 
can be understood using a set of topological Hunds rules. These rules relate the multi-electron 
energy levels to spin and the inter-dot tunneling t, and control charging energies. We map out 
the charging diagram for up to N = 6 electrons and predict a spin-polarized phase for two holes. 
The theoretical charging diagram is compared with the measured charging diagram of the gated 
triple-dot device. 

PACS numbers: 73.21.La,73.23.Hk 
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I. INTRODUCTION 



Following on earlier work which showed that a small and well-controlled number of elec- 
trons can be confined in a single^ and a double quantum dot,—&&-&l an artificial lateral 
quantum molecule consisting of three quantum-mechanically coupled lateral quantum dots 
has been demonstrated.- The triple quantum dot molecule is a natural step toward creat- 
ing quantum dot networks, with potential applications in quantum computing.-^^ When 
filled with three electrons, one electron per dot, this device can serve as a simple quantum 
logic circuit, with each electron spin treated as a qubit. One can also use the molecule as a 
single coded qubit,— ^2.dlM whose states are encoded in the states of three electronic spins 
but tunable with applied voltage. The triple dot could also be used to create entanglement 
between spin qubits,— spin and charge qubits, 16 as a charge rectifier,— 1 ^ and may exhibit 
a characteristic Kondo effect when coupled to the leads.— i22i2L22i22iMi2£i2ii With electrons lo- 
calized on individual dots and their tunneling controlled by gates, the triple dot molecule 
can be also thought of as an implementation of the tunable Hubbard model, an important 
step toward realization of "quantum materials" .2Zi2£i22i22i2ii22 

The electronic properties of the triple quantum dot with one electron per dot have been 
studied theoretically by a number of authors. To make contact with the pairwise-exchange 
formalism used in quantum information,— attempts were made to map the properties of 
this system onto those of the three-spin Heisenberg model. Scarola and Das Sarma^ used 
the Hubbard, variational, and exact diagonalization approaches to demonstrate that this 
mapping can be carried out only for a limited range of triple-dot parameters. Mizel and 
Lidar— arrived at similar conclusions using the Heitler-London and Hund-Miilliken 
schemes to calculate the energy levels of three coupled dots with one electron per dot. 
In both cases the many-body effects were responsible for the appearance of higher-order 
terms in the effective spin Hamiltonian. In an alternative approach, in Ref. [l2| we have 
used real-space wave functions and the configuration-interaction technique to analyze the 
three-electron triple-dot molecule acting as a single coded qubit and shown how its energy 
levels can be tuned by voltages applied to gates defining the structure. 

Properties of the triple-dot molecule as a scattering center have also been studied us- 
ing quantum transport techniques. Using the density-functional and quantum Monte Carlo 
methods, Stopa- 7 calculated the current flowing through a nominally empty molecule con- 



2 



nected to electron reservoirs and under bias. The rectifying behavior of the system predicted 
in this analysis was confirmed experimentally— Landron de Guevara and Orellana^ calcu- 
lated the zero-temperature conductance through a linear molecule coupled in parallel to the 
leads using a Hubbard approach in a magnetic field. Apart from the Fano resonances in 
the spectrum, they found evidence of formation of the quantum-molecular states decoupled 
from the leads. The Hubbard model has also been used to investigate the triple-dot system 
in the Kondo regime, both in the linear— i2L22i and triangular topology— i2i^L2& 

In this paper we describe the electronic properties of a lateral triple quantum dot molecule 
as a function of electron numbers. In analogy to the work on quantum materials,— 1 ^ we 
model our system with the Hubbard Hamiltonian, but the obtained results are verified by 
microscopic methods. In the Hubbard model we retain only one lowest-energy orbital per 
dot. The lowest-energy shell of the molecule can be filled with up to N e = 6 electrons. We 
analyze in detail the ordering of energy levels, the spacing of Coulomb blockade peaks and 
the charging and spin phase diagram of this shell. We demonstrate that the energy levels of 
the molecule are related to the total spin of electrons but not directly related to the charge 
e. We find the spin singlet as the two-electron ground state, with the singlet-triplet (S-T) 
splitting proportional to the single-particle tunneling matrix element t. This is in contrast 
to atoms, where the S-T splitting is proportional to the electronic exchange and hence to e 2 , 
or to magnetic solids, where super-exchange leads to S-T splitting proportional to 1/e 2 . On 
the other hand, for two holes (N e = 4) we predict a spin polarized ground state and a singlet- 
triplet transition driven only by modifying the topology of the system. For three electrons 
in a half-filled shell (N e = 3) we confirm the existence of the frustrated antiferromagnetic 
ground state.— The fact that the tunneling alone distinguishes singlet and triplet states is 
related to the interplay of the Fermi statistics and system topology. We term the set of rules 
established here and relating spin of the ground state to the filling of the shell, topology, 
and tunneling, "topological Hunds rules". The ability to tune tunneling by gates opens the 
possibility of directly manipulating the electron spin using electrical means only, of interest 
in designing novel quantum materials, magneto-electronics and quantum computation. We 
show that the Hubbard model is capable of reproducing the charging diagram of a lateral 
gated triple-dot measured recently by Gaudreau et al. 8 

The paper is organized as follows. In Sec. [Til we describe the model lateral triple-dot 
device and construct the Hubbard Hamiltonian. In Sec. II I II we determine the electronic 
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structure of the device charged with N e = 1 to 6 electrons. Results of the Hubbard model are 
tested against real space (RSP) configuration interaction (RSP-CI) and linear combination 
of atomic or quantum dot orbitals (LCAO-CI) calculations. The charging diagram as a 
function of the dot energies is presented and analyzed in Sec. IIVI In Sec. |V] we relate the 
calculated and measured charging diagrams. Summary and conclusions are presented in 
Sec. ED 



II. THE MODEL 



The proposed model gated triple-dot device realizing the triple dot using only metallic 
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gates, studied in Ref. [12| and related to the one studied by Gaudreau et al. in Ref. 
shown in Fig. [Tj(a). It consists of a heteroj unction with a two-dimensional electron gas 
(2DEG) created at a distance D below the top surface of the sample. The metallic gates 
deposited on the surface serve to deplete the 2DEG underneath. Any opening in the gates 
is translated electrostatically into a local potential minimum, capable of confining a small 
number of electrons. Thus, in our model the three circular holes in the main gate (shown 
in gray) define a triangular triple quantum dot lateral confinement. Each isolated potential 
minimum gives rise to a quantized energy spectrum, of which we retain only the lowest 
energy level E{ in dot i. By tuning the voltage on the main gate we can control the number 
of confined electrons. For example, in Fig. [Jja) we show N e = 2 electrons with parallel spins 
localized on two of the dots. This is not, however, a depiction of a quantum molecular state: 
due to the interdot coupling the electrons are delocalized across the molecule. The main 
gate alone defines a symmetric triangular molecule with identical pairwise coupling of all 
dots. This triple-dot potential can be well approximated by a sum of three Gaussians. 

The single-particle confinement can be additionally tuned by three smaller gates, shown 
in red, green, and blue. Their arrangement with respect to the potential minima is shown 
schematically in Fig. Hfb). The gate VG\ controls simultaneously the lowest energy levels E\ 
and E 2 of dots 1 and 2, and the gate VG 3 controls the energy level E 3 of dot 3. Additionally, 
the gate VGi 3 is designed to tune the topology of the system without significantly changing 
the energies By biasing it with a sufficiently high negative voltage we increase the 
tunneling barrier between dots 1 and 3 and change the sample layout from a closed triangle, 
in which all dots are identically coupled, to a linear molecule, in which the tunneling between 
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dots 1 and 3 is not allowed. 

We examine the electronic properties of our triple quantum dot molecule in the frame 
of the Hubbard model with one spin-degenerate orbital per dot. Without specifying them 
explicitly, the localized orbitals in the Hubbard model are assumed to be orthogonal. This is 
to be contrasted with the approach starting from the linear combination of atomic orbitals 
(LCAO), which are non-orthogonal. The orthogonalization leads to extended, quantum- 
molecular orbitals which serve as a basis for CI calculation. In the Hubbard model, with cf a 
(cj CT ) operators creating (annihilating) electrons with spin o on the orbital of z-th dot, the 
Hamiltonian can be written as: 

3 3 3^3 

H=J2 E i c t c i*+ Z tijcfaCj* + Yl U i n H n H + 2 Z V ueiQj, (!) 

a,i=l a,i,j=l,i^j i=l t,i=l,i^' 

where n ia = cf a Ci a and Qi = riii + n^ are, respectively, the spin and charge density on the z-th 
dot. The above Hamiltonian is characterized by the energy levels of the z-th quantum dot Ei, 
the tunneling matrix elements ty between dots % and j, the on-site Hubbard repulsion Ui, and 
the direct Coulomb matrix elements Vij between dots i and j. These Hubbard parameters 
are schematically shown in Fig. [H(b). With one energy level per dot the triple-dot molecule 
can be filled with up to N e = 6 electrons. 

III. ELECTRONIC STRUCTURE OF THE TRIPLE DOT WITH 1 TO 6 ELEC- 
TRONS 

A. One electron and one hole 

We look for the eigenenergies and eigenvectors of the Hamiltonian (TTJ using the exact 
diagonalization approach. To this end, we create all possible configurations of N e electrons 
on the three localized orbitals, write the Hamiltonian matrix in this basis, and diagonalize it 
numerically. In the simplest case of N e = 1 the basis contains three non-overlapping states, 
{|1), |2), |3)}, where \i) = ct\0) and |0) denotes the vacuum. In this basis the diagonal 
Hamiltonian matrix elements are (i\H\i) = Ei and the off-diagonal elements (i\H\j) = tij. 
With the three dots on resonance, i.e., with E\ = E 2 = E 3 = E and t\2 = £23 — ^13 — t, the 
one-electron energy spectrum is composed of one level with energy Ea — E + 2t, and one 
doubly-degenerate level with energy Eg = Ec = E — t. The order of these levels depends 
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on the sign of the element t. In numerical calculations of the single-particle spectrum 
corresponding to the potential produced by metallic gates 12 shown in Fig. []Ta) we find the 
ground state to be non-degenerate, indicating that t < 0. Additionally, the magnitude of 
the tunneling matrix element can be found from the single-particle energy gap A = 3\t\. 

Knowledge of the sign of the off-diagonal element allows us to construct the single-particle 
molecular orbitals. The ground state is \Mi) — -7= + |2) + |3)), while the two degenerate 
excited states are \M 2 ) = \ (|1) - |2)) and |M 3 ) = (|1) + |2) - 2 • |3». The states \M 2 ) 
and \M 3 ) were chosen to be symmetric with respect to a mirror plane passing through the 
dot 3 and intersecting the (1 — 2) base of the triangle at its midpoint. However, due to the 
degeneracy of the two levels, any pair of orthogonal states created as linear combinations of 
\M 2 ) and \M 3 ) will be viable as eigenstates. The degeneracy of the excited states is a direct 
consequence of the symmetry of the triangular molecule. Changing its topology, e.g., by 
increasing the tunneling barrier between dots 1 and 3, will remove the degeneracy. In the 
limit of an infinite barrier, i.e., £13 = 0, we deal with a linear triple-dot molecule, whose single- 
particle energy spectrum consists of three equally spaced levels: (E — V2\t\, E, E + \/2|£|). 
Thus, the triangular triple dot design makes it possible to engineer the degeneracy of states 
solely by electrostatic means. 

Now we can start to populate our triple-dot molecule with electrons. Let us start our 
many-body analysis with the simplest case of N e = 5. As the maximal number of electrons 
in our system is six, we can interpret the five-electron configurations as those of a single 
hole. The hole (e.g., with spin down) can be placed on either of the dots, and thus our 
basis consists of three configurations: = hf\\N e = 6) = cti c ti c t\ c t\ c i\ |2^) = 

h^\N e = 6) = c^cJcJcJc^lO), and |3 W ) = /iJ|iV e = 6) = c^c J |0), with hf a 

being the creation operator of the hole with spin a on the z-th dot. It is convenient to 
express the energies of these configurations with respect to the total energy of the system 
with six electrons E F = 2E\ + 2E 2 + 2E 3 + Ui + U 2 + U 3 + AV 12 + 4V 13 + AV 23 . We have 
then E[ H) = Ep — E\ — XJ\ — 2V 12 - 2V 13 , E { 2 H) = E F - E 2 - U 2 - 2V 12 - 2V 23 , and 
Ef ' = E F - E 3 - U 3 - 2 Via - 2 V 23 . The three energies are respectively the diagonal terms of 
our single-hole Hamiltonian. The off-diagonal terms are composed out of the single-particle 
tunneling matrix elements. We have (i^ H ^\H\j^ H ^) = — iy-; the negative phase is due to the 
anticommutation relations of the electronic creation and annihilation operators. As we can 
see, the single-hole Hamiltonian can be obtained from the single-electron Hamiltonian by 
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appropriately modifying the diagonal terms and setting Uj <-> —Uj. This is the signature 
of the particle-hole symmetry— However, for the triangular triple dot on resonance this 
symmetry is not reflected in the energy spectrum of the hole: in this case, the opposite sign 
of the off-diagonal element leads to a doubly-degenerate hole ground state. This property is 
immediately apparent in the molecular basis: we create the lowest-energy configuration by 
filling the molecular ground state \M\) with two of the five electrons, and distributing the 
remaining three on the degenerate orbitals \M 2 ) and \M%). The latter can be accomplished 
in two energetically equivalent ways, hence the double degeneracy. Note, however, that the 
electron-hole symmetry is fully restored upon transition to the linear triple-dot molecule. 
For this topology, the single-particle spectrum of both the electron and the hole consists of 
three equally spaced non-degenerate levels. 

B. Two electrons and two holes 

The interplay of topology and statistics is particularly important in the cases of two 
electrons and two holes confined in the triple dot molecule. Let us consider the case of N e = 2 
first. Since the Hamiltonian (JTj) commutes with the total spin operator, we can classify the 
two-electron states into singlets and triplets. Working with the molecular basis set, we form 
the configuration with the lowest energy by placing both carriers with antiparallel spins on 
orbital \M\). Therefore we expect the ground state of the two-electron system to be a spin 
singlet, irrespective of the molecule's topology. However, in order to examine the topological 
and statistical effects in the energy spectrum and the structure of the wave functions, we 
carry out a systematic analysis in the localized basis. 

Due to Fermi statistics, the two electrons with parallel spins cannot occupy the same 
quantum dot. Hence there are only three possible triplet configurations, \Ti) = cJc^jJO), 
\T 2 ) = cJcJjjO), and |T 3 ) = c^cj |0), shown schematically in Fig. [2(a). The three triplet 
configurations interact with each other only via the single-particle tunneling Hamiltonian. 
However, in evaluating the respective matrix elements we need to follow the Fermionic 
anticommutation rules of the creation and annihilation operators. For example, acting with 
H on the configuration T\ to produce the configuration T3 requires the evaluation of the 
following expression: H\Ti) = +£31 Cijejejj |0). In order to remove the electron 1 we first 
have to move it around electron 2, and so H\T\) = — tuct\ c ti c il c t\ |0) = — ^31 IX3) - Hence, 
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tunneling of the electron from dot 1 to dot 3 in the presence of the electron in dot 2 generates 
an additional phase or changes the sign of the tunneling matrix element. This is of course 
the most elementary property of Fermions brought out so clearly in this simple model. By 
contrast, tunneling from dot 2 to dot 3 in the presence of electron in dot 1 does not change 
the sign of the tunneling matrix element. The resulting triplet Hamiltonian matrix takes 
the following form: 



H' 



T 



El + E% + V\2 ^23 — ^13 

^23 E\ + £?3 + V13 t\2 

— *13 tl2 E 2 + E 3 + V 2 3 



(2) 



Ht is related to the one-hole Hamiltonian. This similarity becomes more apparent if Ht is 
written in the basis { | T\ ) , — |T 2 ), |T 3 )}, in which case all the off-diagonal elements acquire a 
negative phase. This is not surprising, since the single- hole configurations analyzed in the 
previous Section can be generated from the above triplet configurations simply by adding 
to them an inert core of three electrons spin up, one electron per dot. With the three dots 
on resonance and all tunneling matrix elements tij equal and negative, the triplet energy 
spectrum is found to be (2E + V — \t\, 2E + V — \t\, 2E + V + 2\t\). As in the case of the single 
hole, the lowest-energy triplet state is doubly degenerate. Moreover, the renormalization of 
the lowest energy 2E + V — \t\ from the single configuration energy 2E + V, as well as the 
gap in the triplet spectrum, are determined entirely by tunneling. The splitting between 
the ground and first excited states is the same as that found in the single-carrier case and 
equals 3\t\. 

We shall now demonstrate that topology and statistics differentiates between triplet and 
singlet two-electron states. The singly-occupied singlet configurations \Si), |S 2 ), and IS3) 
are obtained from the triplet configurations |Ti), |T 2 ), and |T 3 ) by flipping the spin of one 
electron and properly antisymmetrizing the configurations. For example, the configuration 
I £1) = -7= (c^c^ + c^c^) |0). In addition to the singly-occupied configurations there are 
also three doubly-occupied configurations, e.g., |S 4 ) = c]jC^|0), as shown in Fig. 12(b). In 
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the basis of the six configurations the two-electron singlet Hamiltonian can be written as: 



Ei + E 2 + V12 


^23 


tl3 


V2t 12 







^23 


E± + E 3 + V13 


tl2 


y/2t 13 







tl3 




E 2 + E 3 + V 23 





a/2*23 




V2t 12 







2E 1 + C/i 








V2t 12 





V2t 23 





2£ 2 + C/ 2 










V2t 23 








2£ 3 + U 3 



• (3) 



The 3x3 upper left-hand corner of H s corresponds to the three singly occupied configura- 
tions I Si), |S 2 ), and |S 3 ). It is similar to the two-electron triplet Hamiltonian H T but differs 
from it by the positive phase of the tunneling matrix element t 13 . Hence, in the triangular 
topology of the triple-dot molecule the tunneling from dot 1 to dot 3 distinguishes between 
the singlet and the triplet spin configurations. By setting t± 3 = 0, i.e., upon transition to 
the linear topology, this difference disappears. However, the singlet basis is still different 
from its triplet counterpart due to the presence of the doubly-occupied configurations. 

For the dots on resonance the energies of the six singlet levels can be obtained analytically. 
The spectrum can be grouped into two non-degenerate levels Ef 2 : 



Ef 2 = (2E + V -2\t\) + - 



{U-V+ 2\t\) ± ^(AV2ty + (y -v + 2\t\y 

and two groups of doubly degenerate levels Ef_ 6 : 



(4) 



El_ G = (2E + V + |i|) + \ (U-V- \t\) ± y/(2y/2t) 2 + (U-V- \t\) 2 



(5) 



In the strong coupling limit U ^> V > \t\ the singlet ground-state energy Ef ~ (2E + 



V - 2\t\) - 



8t 2 



r _ v -, while the triplet energy Ef = (2E + V — \t\). Thus, the two-electron 
ground state is always a spin singlet. The singlet-triplet gap, separating Ef from E\ is 



t + 



St 2 
U-V 



It is proportional to the tunneling matrix element \t\ and contains the 



second-order super-exchange correction ~ t 2 / (U — V) due to the doubly occupied singlet 
configurations. Removing the resonance by detuning the onsite energies E^ enhances the 
contribution from the doubly-occupied states. Therefore the ground state maintains its 
singlet character independently of the choice of gate voltages. 

The situation is qualitatively different when two holes, instead of two electrons, populate 
the system. The two holes are created when two electrons are removed from the closed-shell 
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configuration with N e = 6, i.e., they correspond to N e = 4 electrons. In the molecular basis 
corresponding to the triangular triple dot we put two electrons on the lowest-energy orbital 
| Mi), and the remaining two electrons on the degenerate pair of orbitals |M 2 ) and |M 3 ). 
With this alignment of levels it is possible to create both triplet and singlet configurations, 
all with the same single-particle energy, and it is not immediately clear which total spin is 
preferred. On the other hand, in the limit of the linear triple dot the molecular orbitals are 
non-degenerate and the four-electron ground state is expected to be a spin singlet. 

The selected two-hole singlet and triplet configurations in the localized basis are illus- 
trated in Fig. E(c). Let us focus on the triplets first. They involve one electron spin-up 
occupying the first, second, or third dot in the presence of an inert core of three spin- 
down electrons. For example, the configuration shown in left-hand panel of Fig. [2](c) 
can be written as {T^) = hf\h2\\N e = 6) = ci-CgicJ cf, |0). Therefore, the hole triplet 
Hamiltonian is equivalent to the single-electron Hamiltonian, differing from it only in di- 
agonal terms. For example, the energy of the configuration {T^) is (Ti\H\Ti ) = 
Ep — Ex — E2 — U\ — U2 — 3Vi2 — 2V13 — 2V23. The two-hole triplet Hamiltonian can 
also be compared to the two-electron triplet Hamiltonian Ht, written in the modified ba- 
sis set {\Ti) , —\T 2 ) , \T 3 )} (i.e., with all off-diagonal matrix elements acquiring a negative 
phase). Setting aside the diagonal matrix elements, the two Hamiltonians are connected 
by the electron-hole symmetry transition t^- <-> However, unlike that of the elec- 

tronic triplet, the ground state of the hole triplet is non-degenerate, and its energy is 
Ex = Ep — 2E — 2U — 7V — 2\t\. As it is in the case of the single electron and the 
single hole, the particle- hole symmetry between the two-electron triplet and the two- hole 
triplet is fully restored upon transition to the linear topology of the triple dot. 

Let us move on to considering the two-hole singlet configurations. The singly-occupied 
states, illustrated in the middle panel of Fig. M^c), involve the two holes occupying two 
different dots, while the doubly-occupied states, such as the one in the right-hand panel of 
Fig. [2(c), hold both holes on the same dot. The two-hole singlet Hamiltonian is analogous 
to that of the two-electron singlet, Eq. (J3]). However, we need to replace the energy of 
two-electron complexes with the energy of two-hole complexes, and change the phase of the 
off-diagonal elements connecting the singly-occupied configurations. The sign of elements 
connecting the singly- and doubly-occupied configurations does not change, which 
breaks the particle-hole symmetry. 
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The ground-state energy of the hole singlet for the triangular triple dot on resonance 
is well approximated by E? (H) « (E F - 2E - 2U - IV - \t\) - Compared to the 

energy of the triplet, E 1 is increased by the tunneling element \t\, but decreased by 
the super-exchange contribution jj^y- Note that the two-hole super-exchange term is four 
times smaller than the super-exchange correction to the energy of the two-electron singlet. 
By increasing the tunneling or decreasing the on-site Hubbard repulsion we can increase 
the contribution from super-exchange and lower the energy of the singlet state. Therefore, 
the total spin of the two-hole ground state for the triple dot on resonance depends on the 
interplay of Hubbard parameters. For 2\t\ < U — V the ground state is a spin triplet, 
and a triplet-singlet transition can be induced by increasing the hopping matrix element. 
The triplet-singlet transition can also be induced by biasing one of the dots, which lowers 
the energy of the doubly-occupied singlet configurations. Hence the configuration of two 
holes shows a nontrivial dependence on tunneling, Coulomb interactions and gate voltages 
allowing to control the system's magnetic moment purely by electrical means. 

C. Three electrons 

To complete our understanding of the energy levels of a triple quantum dot molecule we 
need to analyze the half-filled case of three electrons (or, equivalently, three holes). We start 
with the completely spin-polarized system, i.e., one with total spin 5 = 3/2. In this case 
we can distribute the electrons on the three dots in only one way: one electron on each site, 
which gives a spin-polarized state |a 3 / 2 ) = c ^ c i c ii|0}- As the basis of our Hilbert space 
consists of one configuration only, |a 3 / 2 ) is the eigenstate of our system, and its energy is 
£3/2 = Ei + Ei + E 3 + V12 + V13 + V23. Let us now flip the spin of one of the electrons. 
This electron can be placed on any orbital, and with each specific placement the remaining 
two spin-down electrons can be distributed in three ways. For example, Fig. [31(a) shows the 
three configurations with the spin-up electron occupying the dot 1. Thus, altogether we can 
generate nine different configurations. Three of these configurations involve single occupancy 
of the orbitals. They can be written as \a) = c^c^cfAO), \b) = ct\ c t\ c 2-\\^) ■> anc ^ l c ) = 
cjc^i C3J0}. Out of these three configurations we construct the three eigenstates of the total 
spin operator. One of those eigenstates is |a 3 / 2 ) = ^(\ a ) + \b) + \c)) , and it corresponds to the 
total spin 5 = 3/2. The two other eigenstates, \ai/ 2 ) = -4^(|a) — \b)) and | &1/2) = 77g(l a ) + |fr) — 
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2 1 c) ) , correspond to the total spin S = 1/2. The remaining six configurations involve doubly- 
occupied orbitals. They are |c 1/2 ) = cjcjc^|0), \d 1/2 ) = cjcjcj |0), |ei /2 ) = cJc^cJ T |0), 

I/1/2) = c u c 2j. c 2tl0)> Idi/z) = c U c 3i c 3T |0) ' l^i/a) = c i c 3i c 3t|0)- All these configurations are 
eigenstates of the total spin with S = 1/2. Thus, among our nine spin-unpolarized states 
we have one high-spin, and eight low-spin states. In this basis the Hamiltonian matrix is 
block-diagonal, with the high-spin state completely decoupled. The energy corresponding 
to this state is equal to that of the fully polarized system discussed above, and is equal to 
-E3/2. In the basis of the nine S = 1/2 configurations we construct the Hamiltonian matrix 
by dividing 9 configurations into three groups, each containing one of the singly-occupied 
configurations \a), \b), and |c). By labeling each group with the index of the spin-up electron, 
the Hamiltonian takes the form of a 3 x 3 matrix: 

#1/2 = 

The diagonal matrix, e.g., 

2E 1 + E 2 + 2V 12 + Ux ^23 ~~ ^13 

t 23 2E 1 + E 3 + 2V 13 + f/i to 

— £13 £12 Ei + E 2 + E 3 + Vi 2 + V%3 + V 23 

describes the interaction of three configurations which contain spin-up electron on site 1, i.e., 
two doubly-occupied configurations lci/2) and \di/ 2 ), and a singly-occupied configuration \a) 
(in this order, see Fig. [3](a)). The configurations with double occupancy acquire the diagonal 
interaction term U. The three configurations involve a pair of spin-polarized electrons (spin 
triplet) moving on a triangular plaquette in the presence of a "spectator" spin-up electron. 
Because of the triplet character of the two electrons, the phase of the hopping matrix element 
— £13 from site 1 to site 3 is different from the phase of the hopping matrix element -K23 from 
site 2 to site 3. As discussed above, the negative phase in — 1% 3 distinguishes the singlet and 
triplet electron pairs. The remaining matrices corresponding to spin-up electrons localized 
on sites 2 and 3 can be constructed in a similar fashion. The interaction between them is 



#1 


f l2 


f+ 
J 31 


J 12 


H 2 


^23 


T 3 1 


f+ 


H 3 



(6) 
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given in terms of effective hopping matrix 



T 



o -u 





+tij 










-t 





1.1 



There is no direct interaction between the configurations with single occupancy, since such 
scattering process would have to involve two electrons, one with spin up and one with 
spin down. This cannot be accomplished by the single-particle tunneling. These states are 
coupled only indirectly, involving the configurations with double occupancy. 

The low-energy spectrum of the Hubbard Hamiltonian of the N e = 3 quantum-dot 
molecule can be further approximated by the spectrum of the model spin Hamiltonian : 



H 3e = E 3/2 + J l3 (Si ■ S, - 1/4) + D ^Si ■ (S 3 x S k ) 

i<j<k 



(7) 



Here, E 3 / 2 is the energy of the spin S — 3/2 stcitc, Jij are exchange matrix elements of the 
Heisenberg part of the spin Hamiltonian which depend on microscopic parameters of the 
triple dot, and Djjk are higher order spin-spin interactions discussed, e.g., by Scarola and 



Das Sarma in Ref. 
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We define the effective exchange constant J for the triple dot molecule with three electrons 
in terms of the gap between the S = 1/2 and 5 = 3/2 states as £3/2 — Eu 2 = 3 J/2. Without 
the higher order corrections J would have been equal to the Heisenberg J^, otherwise it is 
simply related to the gap of the N e = 3 electron spectrum. 

The mapping of the behavior of our system onto the effective exchange Hamiltonian 
(J7J) connects our analysis to the general formalism used in quantum computing.— 1 ^ Our 
considerations do not introduce any new elements into that formalism, but rather provide 
means for its realistic and accurate parametrization, reflecting the properties of an actual 
gated triple-dot device. 

The Heisenberg Hamiltonian ([7j) can be used to model the behavior of three electrons 
confined in a triple dot treated as three coupled qubits. However, it applies also to a coded 
qubit scheme, in which the states of the entire molecule are treated as the logical states of a 



single qubit. In Ref. [l2| we have presented a detailed analysis of such a system, in which we 
selected the two lowest total spin 1/2 states as the logical states |0i) and |1^) of the coded 
qubit, respectively. In that design, the control of the energy gap between the two states by 
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the gate voltage provides means for the single-qubit operations. Again, our current work 
allows to parametrize this model with the Hubbard parameters appropriate for a specific 
triple quantum dot. 

D. Comparison of Hubbard, LCAO-CI and RSP-CI results 

We shall now find the values of the Hubbard parameters appropriate for a typical triple 
quantum dot system. These parameters are obtained by fitting the electronic properties 
discussed above either to results of microscopic calculations, or to experimental data. In 
this Section we will focus on the former, while the latter will be discussed in Section |V] 

In what follows we shall express all energies in units of the effective Rydberg, 171 = 
m* e 4 / '2e 2 h 2 , and all distances in units of the effective Bohr radius, lag = eh 2 /m*e 2 , where 
e and m* are the electronic charge and effective mass, respectively, and e is the dielectric 
constant of the material. For GaAs parameters, m* = 0.067mo and e = 12.4, we have 
17?. = 5.93 meV and lag = 9.79 nm. As the model lateral triple-dot system we take the 
structure shown in Fig. Ufa), discussed by us in detail elsewhere.— We take the main gray 
gate to be a square with the side length of 22. 4ae- The diameter of each circular opening 
is 4.2a#, the distance between the centers of each pair of the holes is 4.85ob- The gate is 
positioned 14ob above the two-dimensional electron gas and a voltage of — \e\V = 107?. is 
applied to it to create the symmetric triangular triple quantum dot. 

We focus on the case of N e = 3 confined electrons. Our analysis consists of two steps. 
First, we find N$ lowest-lying single-particle energies and wave functions of the system and 
obtain the Coulomb matrix elements involving all these states. Second, we calculate the 
three-electron eigenenergies within the configuration-interaction (CI) approach. 

The one-electron properties of the system can be derived in a real-space approach (RSP) 
involving numerical diagonalization of the discretized single-particle Hamiltonian.~ We com- 
pute N s = 9 lowest-lying levels. The ground state is separated from the first excited state 
by an energy gap of 0.18777?., while the gap between the first and second excited states is 
much smaller and equal to 0.00617?.. This agrees well with the Hubbard model, predicting 
a degeneracy of the two excited states. Also, from the average gap between the ground 
and excited states, which in the Hubbard model equals 3\t\, we can extract the tunneling 
parameter tyi = £13 = t 2 3 = —0.06367?.. 
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The RSP approach, while being accurate, is computationally intensive. As an alternative 
we consider a method based on the linear combination of atomic orbitals (LCAO).— To this 
end, we approximate the numerical triple-dot lateral confinement, obtained as a solution of 
the Poisson equation, with a sum of three Gaussians: 



The pairs (xi,yi) are coordinates of the center of each dot. For our symmetric triple dot 
a good fit is obtained for d x = d 2 = d 3 = 2.324a B and V {1) = V {2) = v} 3) = 5.86411. 
We seek the quantum-molecular single-particle states in the form of linear combinations of 
single-dot orbitals localized on each dot. To simplify the calculations, we take these orbitals 
to be harmonic-oscillator (HO) wave functions of a two-dimensional parabolic potential, 
obtained by extracting the second-order component from each Gaussian. We take one s- 
type HO orbital per dot, and solve the generalized single-particle eigenproblem formulated 
in this nonorthogonal basis set.— As a result, we obtain N$ = 3 quantum- molecular levels: 
a non-degenerate ground state and a doubly degenerate excited state, separated by a gap of 
0.03547?.. This structure of levels is reproduced by the Hubbard model with the tunneling 
parameter t = —0.01187?. Note that in the LCAO case the tunneling gap is much smaller 
than that obtained in the RSP calculation. This is due to the restricted LCAO basis set, 
which underestimates the overlap between the single-dot orbitals. The agreement between 
the two approaches can be improved upon inclusion of the p and d HO orbitals in the LCAO 
basis, at the expense of clarity. 38 

With the single-particle energies and the Coulomb matrix elements calculated using the 
quantum-molecular orbitals, we can now proceed to the CI calculation of three-electron 
properties. We create all possible configurations of the three electrons with total S z = —1/2 
on Ns quantum-molecular states (Ns = 9 for the RSP, and Ns = 3 for the LCAO approach), 
build the many-body Hamiltonian matrix in the basis of these configurations, and diagonalize 
it numerically— The resulting spectra are shown in the left-hand parts of Fig. [3](b) for the 
RSP-CI, and in Fig[3^c) for the LCAO-CI approach. From the gaps separating the levels we 
can extract the Hubbard interaction parameters. They are: V = 0.4797? and U = 1.5397? 
in the RSP-CI case, and V = 0.4227? and U = 2.55 77? in the LCAO-CI case. The resulting 
Hubbard three-electron spectra are plotted in the right-hand parts of the Figures [3](b) and 





(c). 
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All models predict the ground state of the system to be doubly degenerate and to have 
total spin S = 1/2. This is easily understood by building the lowest-energy configuration 
with triple-dot molecular orbitals: two out of three electrons are placed on the orbital \M\) 
with antiparallel spins, and the third electron - on one of the degenerate orbitals \M 2 ) or 
|M 3 ). This configuration has total spin S = 1/2, and no spin transition is expected upon 
the change of the system's topology. Further, in all spectra the first excited state has total 
spin S = 3/2. This is the (03/2) state from the previous Section, equivalent to the spin- 
polarized configuration with one electron per dot. The energy gap between the low-spin and 
the high-spin states, expressed in the language of the effective Heisenberg Hamiltonian, is 
equal to 3 J/2, with J > (the ground state is antiferromagnetic). 

The remaining excited states, involving doubly occupied configurations, are visible at 
higher energies. They are separated from the low-energy, singly-occupied states by a gap 
proportional to the Hubbard onsite interaction parameter U . This parameter is larger in 
the LCAO-CI approach because of the relatively small spatial extent of the HO basis states. 
This is consistent with the underestimated tunneling gap found earlier in the calculation of 
the single-particle spectra. On the other hand, the Hubbard interdot interaction parameter 
V is similar in both approaches. In the Hubbard model, the high-energy part of the spectrum 
is composed of three doubly-degenerate states. The degeneracy of the lowest and the highest 
level within this band is well reproduced in both microscopic models, while the middle level 
appears to be split by a smaller gap in the LCAO-CI, and a larger gap in the RSP-CI 
spectrum. Finally, the RSP-CI result reveals further levels, with both total spin 5=1/2 
and S = 3/2. Their appearance is a consequence of the extended basis, containing Ns = 9 
single-particle molecular states, compared to Ns = 3 states in LCAO-CI and Hubbard 
models. The additional states can be accounted for systematically by including more than 
one orbital per dot in the localized basis set.— To conclude this analysis, we find the Hubbard 
model to give qualitatively correct results but caution has to be exercised when making a 
quantitative comparison. 

Up to now we have explored the case of a symmetric triangular triple-dot molecule. Let 
us now tune the topology of the system using the gate GV13. Figure 0(a) shows the three- 
electron spectrum obtained with the RSP-CI method as a function of the voltage applied to 
this gate, and the low-energy part of this spectrum is shown in the inset to this Figure.— 
In the Hubbard model, this change of topology can be accounted for by tuning the single- 
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particle tunneling parameter £13. The corresponding spectra are shown in Fig. H](b). Both 
approaches indicate that the change of topology of the molecule leads to a splitting of the 
two degenerate 5 = 1/2 levels. This property is a consequence of the removal of degeneracy 
of the single-particle molecular orbitals: with two electrons forming a spin singlet on the 
orbital |Mi), the third electron probes the splitting between the levels |M 2 ) and |M 3 ) The 
ability to tune the splitting by electrostatic means only suggests a possible use of the two 
low-spin states as logical states of a voltage-controlled coded qubit.— In the language of 
Hubbard configurations discussed in the previous Section, these states can be written as 
|0 L ) = ao^(|a) - +/3 |A ) and \1 L ) = «i^(|a) + |6) - 2|c>) +/3 a | A x ), where |A ), |A X ) 
are contributions of the doubly-occupied configurations. 



IV. CHARGING DIAGRAM OF THE TRIPLE DOT 

We can now construct the charging diagram of the triple-dot molecule. For any num- 
ber of electrons N e (1 to 6) and any pair of gate voltages, or equivalently, quantum dot 
energies E i: we can establish the ground-state energy E GS (N e ) by diagonalizing the Hub- 
bard Hamiltonian. We use these energies to calculate the chemical potential of the triple 
quantum dot molecule fJ,(N e ) = E GS (N e + 1) — E GS (N e ). When (J,(N e ) equals the chemical 
potential fi^ of the leads, the N e + 1st electron is added to the iV e -electron quantum-dot 
molecule. This establishes the total number of electrons N e in the quantum dot molecule 
and their total spin as a function of applied voltages, or quantum-dot energies. Changes in 
electron numbers can be detected by the Coulomb blockade (CB), spin blockade, or charg- 
ing spectroscopies.-^ The calculated stability diagram, with Hubbard parameters extracted 
from the RSP-CI calculation for three electrons, is shown in Fig. [5j 

Figure [5(a) shows the addition spectrum for the triple dot on resonance, i.e., when all 
dots are characterized by the same onsite energies, tunneling amplitudes, and Coulomb 
matrix elements. We follow the addition spectrum as we change the onsite energy E of each 
quantum dot with respect to the chemical potential of the leads \xl = 0. From the condition 
fiL = E — 2\t\, the energy E(l) corresponding to the addition of the first electron equals 
twice the hopping matrix element, E(l) = 2\t\. At this energy the first Coulomb blockade 
peak of the triple quantum dot molecule should be observed. 

The onsite energy corresponding to the second CB peak, i.e., when the second electron 
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enters the dot, equals E{2) = —V + jjzy- The energy to add the second electron, or 
the spacing between the first two peaks A 12 — V + 2\t\ — -jjzy, is proportional to the 
direct Coulomb interaction V between two electrons on two different dots and to twice the 
tunneling matrix element, and is reduced by super-exchange interaction. 

The third electron enters the molecule for E(3) = —2V + ^ — 2\t\ — and the spacing 
of the third and second CB peak equals A 23 = V + 2\t\ — ^- + 2-^y. This spacing is 
proportional to V, 2t, and twice the super-exchange, but is reduced by the spin gap of the 
N e = 3 electron complex, equal to 3J/2. The difference between the spacing of the (2,1) 
and (3,2) peaks, A 2 3 — Ai 2 = — ^ + directly measures the difference between the 

exchange in the triply occupied quantum dot molecule and three times the super-exchange 
in a doubly occupied quantum dot molecule. 

The half-filled molecule can also be probed by adding the fourth electron. This electron 
enters the dot for E(A) = -U - 2V - ^ + 2\t\. Since the four-electron states are the first to 
be built by doubly-occupied configurations, the corresponding CB peak is spaced from the 
one for the third electron by a large on-site Coulomb energy U. The separation between CB 
peaks equals A 43 = U + 2^f — 4|£| — -fjzy- It reflects the triplet state of two holes, and is a 
measure of U, J, and \t\ but not directly V. 

The expressions for peak spacings can be used to extract the Hubbard parameters of 
the system from the measured CB spectrum of the triple dot on resonance. With these 
parameters we can now explore the full charging diagram of the molecule as a function of 
the onsite dot energies. If each of the energies Ei can be varied independently, the resulting 
stability diagram is three-dimensional, and therefore difficult to visualize. This is why in the 
proposed device, shown in Fig. [U^a), the dots one and two are tuned by a single gate VG\ 
while dot three is tuned by gate VG 3 . Figure [5](b) shows the corresponding cross-section 
of the stability diagram, calculated with the RSP-CI Hubbard parameters. The diagram 
shows the regions (E\ = E2, E%) where different electron numbers are stable. The regions 
are denoted by (Aq, N 2 , iV 3 ) where iVj is the number of electrons (for N e < 3) or holes (for 
N e > 3) in the i-th dot. For example, (1, 1, 1) denotes the half-filled triple dot with one 
electron in each dot. Additionally, the regions are color-coded to indicate the total spin of 
the molecule. We find that the two electrons always form a spin singlet, but the total spin of 
the two- hole system can be changed from the triplet, which is stable close to the resonance 
condition, to a singlet. This transition can be induced by tuning the gate voltages and does 
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not require the presence of a magnetic field. 



V. COMPARISON OF THEORETICAL AND EXPERIMENTAL CHARGING DI- 
AGRAMS 

We now turn to the comparison of theory with experiment. The addition diagram of 
a triple-dot lateral device, measured recently by our group,- is shown in Fig. [S^a). The 
layout of the metallic gates composing the device and the resulting potential minima are 
shown in the inset to this Figure. A similar arrangement of five bottom gates and one 
top gate has been used to define electrostatically the lateral double quantum dot, with the 
quantum point contact (QPC) used as a charge detector.-^ To this end, a sufficiently large 
negative voltage is applied to the gates IB and 5B, as well as to the top gate (T) and the 
middle bottom gate (the gate 3B). The smaller bottom gates, 2B and 4B, are then used as 
plungers, i.e., tuned with a smaller voltage to influence each of the dots locally. However, 
if a large negative voltage is applied to the small gates 2B and 4B, and a smaller voltage 
- to the middle gate 3B, a structure of three potential minima is created: the dot 1 close 
to the gate IB, the dot 3 - near the gate 5B, and the dot 2 - in the middle, between the 
gates T and 3B. In this arrangement, the dots form a linear chain, so the electrons cannot 
tunnel from the surrounding 2DEG directly to the middle dot. In the experimental addition 
spectrum, however, three sets of lines with distinct slopes are detected, indicating that each 
dot is connected to the leads independently. Moreover, two sets of lines show a stronger 
dependence on the voltage V\b than on V 5 b, while one set exhibits an opposite tendency. 
This suggests a formation of a ring- type arrangement, consisting of two dots contained on 
the left-hand side, and one - on the right-hand side of the device, as shown schematically in 
the inset to Fig. El(a). The double potential minimum on the left-hand side is created most 
likely by a mesoscopic fluctuation of the background potential of the sample. This makes it 
difficult to control the dots 1 and 2 independently. On the other hand, the proposed sample 
layout, shown in Fig. [T^a), results in the formation of electrostatically defined triangular 
triple-dot confinement. The design has been adapted to approximate the functionality of 
the experimental device, but it can be modified to allow for independent control of both 
onsite and tunneling energies, however at the expense of a more difficult to fabricate, vertical 
structure of multiple gates. This is why we do not suggest this sample layout as a practical 
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gating scheme, but use it to demonstrate the degree of control we aim to achieve in our 
future designs of triple-dot systems. 

In our experiment, the energy landscape of the lateral confinement is tuned by all six 
gates, but the addition diagram is measured only as a function of two gate voltages, V\b 
and V$b- When these voltages are set to large negative values, the system is completely 
depleted of electrons. This corresponds to the bottom left-hand region (0, 0, 0) of the stability 
diagram. In the Hubbard model this region would correspond to the onsite dot energies Ei 
being larger than the chemical potential of the leads, i.e., the upper right-hand corner (0, 0, 0) 
of the diagram in Fig. E^b). As the gate voltages are made less negative, the energies Ei are 
lowered, and subsequent electrons enter the molecule. These addition events are detected 
in the QPC current Iqpc, reacting to the changes in the charge distribution in the system. 
The dark lines in Fig. HU(a) denote the boundaries between regions corresponding to different 
stable electron numbers. Let us focus on the addition line composed of sections A, B, and C, 
which marks the addition of the first electron to dot 1, 2, or 3, respectively. The quantum 
molecular character of the system is revealed by the curvature of this line close to the regions 
denoted as D and E, where the dots 1 and 2, or 2 and 3 are on resonance, respectively. It 
is clear that the dots 1 and 2 are coupled much more strongly than the dots 2 and 3. The 
dashed lines drawn in these regions connect the points of inflection of the addition lines 
of the first and second electron. In the absence of quantum tunneling, these lines would 
correspond to the conditions E\ = Ei and Ei = E$, respectively. Away from these regions 
the dots are far from resonance, and the electrons are added to orbitals well-localized on 
individual dots. Therefore, the asymptotes drawn with respect to the sections A, B, and C 
of the one-electron line will define the respective single-dot properties. 

As a first approximation the onsite energies Ei are expressed as linear functions of the 
two gate voltages: 

E i = a i V 1 B + PiV B B + 'yi. (9) 

Let us first focus on establishing the coefficients a« and In general, we seek six co- 
efficients, but have only five equations at our disposal (the asymptotes to sections A, B, 
and C, and the two resonance conditions D and E), so at least one parameter has to be 
established independently. In this case, however, from independent measurements we know 
three coefficients: «i = a 2 = —100 meV/V, and /3 3 = —100 meV/V. The equality of the 
coefficients a± and ai is reflected by the vertical character of the line D. The remaining three 
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coefficients can now be easily extracted from the asymptotes to the sections A, B, and C, 
and are ft = -19.0 meV/V, ft = -26.923 meV/V, and a 3 = -22.923 meV/V. 

We are now in a position to convert the charging diagram from the coordinates (Vib, V^b) 
to the energy coordinates, assuming for the moment that the coefficients 7, are zero. Since 
the Hubbard parameters are of one-dot and two-dot character only, they can be systemati- 
cally fitted by extracting the features involving the dot pair from the charging diagram 
and replotting them as a function of the energies (Ei,Ej). As an example we discuss the 
case of % — l,j = 2, with the translated charging diagram shown in Fig. E^b). In this 
Figure, the solid black lines show the experimental data, while the dotted horizontal and 
vertical lines are the asymptotes. As already discussed, in the upper right-hand region, 
corresponding to large values of E\ and E 2 , the system is empty. Starting in this region, 
we can decrease the energy E\ while maintaining E 2 constant: this corresponds to moving 
horizontally across the diagram. Along the way we shall first cross the rightmost vertical 
asymptote, which will mark the addition of the first electron to dot 1 in the zero-coupling 
regime (i.e., t\ 2 = 0), thereby driving the system into the configuration (1,0,0). Because 
of our assumption of the chemical potential of the leads /il = 0, the energy E±, which this 
asymptote defines, is equal simply to —71. In the similar fashion, from the region (0,0,0) 
we can move vertically downwards, decreasing E 2 while keeping E\ constant. Crossing of 
the top horizontal asymptote marks the addition of the first electron into the second dot, 
i.e., formation of a configuration (0, 1, 0), and defines the parameter 72. The two asymptotes 
cross at a right angle, which would be an expected behavior of the addition lines at zero 
coupling. However, the experimental data trace a hyperbola, whose curvature is a direct 
measure of the single-particle tunneling element t 12 . 

Let us now position ourselves in the region in which the first electron has entered the dot 
1 (the region (1,0,0), the top part of the diagram). As we move vertically downwards, we 
encounter the top horizontal asymptote. This line would mark the addition of the second 
electron, and its placement on the second dot, but only in the case the electrons were not 
interacting. Since it is necessary to compensate for the Coulomb off-site charging energy, 
the actual addition takes place at lower energy E 2 , i.e., upon crossing of the horizontal 
asymptote second from the top. The energy distance between the two horizontal asymptotes 
corresponds directly to the Hubbard parameter Vi 2 . An identical value is obtained by 
performing an analogous analysis starting in the region (0, 1, 0) (the right-hand edge of the 
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diagram), and moving horizontally to the left. Finally, we can find the onsite Coulomb terms 
U\ and U 2 by examining the energy differences between asymptotes marking the addition of 
the second and the third electrons. These terms define the size of the stability region (1, 1, 0), 
as shown in Fig. E^b). By using a similar analysis for the features involving the second and 
third dots, and then the first and third dots, we can systematically extract all Hubbard 
parameters. In our case they are (in meV): 71 = —34.238, 72 = —37.169, 73 = —36.246, 
t 12 = -0.053, £13 = t 23 = -0.0077, V 12 = 0.4623, V 13 = 0.0448, V 23 = 0.0962, U x = 2.238, 
U 2 = 2.1262, and U 3 = 1.8923. Figure [TJ^a) shows the charging diagram computed with the 
Hubbard model with the above parameters as a function of the gate voltages. It coincides 
exactly with the experimental diagram shown in Fig. [6](a). 

Note that with the dependence of the onsite energies E^ on gate voltages defined in 
Eq. Q), the triple-dot molecule is on resonance only for (Vi, V5) = (— 0.363U, — 0.37U). In 
Fig. [6](a) this point is found in the (0, 0, 0) region, and this is why the charging diagram 
is essentially a superposition of two double-dot diagrams, and no features unique to the 
resonant triple-dot molecule are visible. It has been demonstrated 8 that by retuning the 
gates making up the device the point of resonance can be shifted to the region of the 
diagram where the electrons start populating the system. This results in the appearance of 
the quadruple points, in which four different electronic configurations are on resonance, and 
charge redistribution effects similar to those in quantum cellular automata. 

Now let us assume that we can control the three quantum dot energies in our device 
independently, with the dots one and two on resonance. Using the Hubbard parameters 
found for our device we compute the charging diagram as a function of E\ = E 2 ,E 3 , and 
shown it in Fig. [T)^b) . The computed charging diagram is similar to that in Fig. [5](b), in 
which we resolve the spins of the electronic states. The results agree with our theoretical 
predictions, including the existence of the triplet four-electron phase. We find this phase 
stable across only a very small range of onsite energies. Detecting the electrostatically driven 
triplet-singlet transition will be investigated in the future. 

VI. CONCLUSIONS 

In conclusion, we presented a theory of electronic properties of a triple quantum dot 
molecule. The electronic properties can be understood in terms of a topological Hunds rule, 
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which determines the spin of the molecule as a function of the filling of the electronic shell. 

When the three dots form a symmetric triangular molecule on resonance, the ground 
state for two electrons is a spin singlet, for three electrons (half-filled shell) it is an anti- 
ferromagnetic 5 = 1/2 configuration, and for two holes it is a triplet. The topology and 
statistics enter through the dependence of the energies of states on total spin. For example, 
the singlet-triplet splitting is found to depend on tunneling and not on charge. The energetics 
and the charging diagram are mapped out, compared with experiment and analyzed in detail. 

We have also demonstrated that the Hubbard model is capable of reproducing the experi- 
mental addition spectra in a quantitative manner. We have described a systematic procedure 
of extracting the Hubbard parameters from the elements of the measured charging diagram. 
Since in the experiment the single-particle orbital energies are controlled by gate voltages, it 
should be possible to induce the triplet-singlet transition for a four-electron molecule purely 
by electrostatic means. Our preliminary calculations indicate that such a transition should 
be possible in the case of the experimental lateral triple-dot device used by Gaudreau et al— 

Future work, including requirements for new device layout, is outlined. 
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FIG. 1: (Color online) (a) Cross-sectional view of a model of the three coupled gated lateral 
quantum dots. The grey rectangular gate contains three circular openings, which translate into 
minima of the electrostatic potential at the level of the two-dimensional electron gas. The red and 
green gates can be used to shift the potential minima of the dots underneath them with respect to 
the rest of the system. The blue gate is used to tune the tunneling barrier between dots 1 and 3. 
(b) Schematic representation of the triple dot structure. 

FIG. 2: (Color online) The three triplet configurations in a two-electron triple dot molecule (a) 
and examples of singly and doubly occupied singlet configurations for two electrons (b) and two 
holes (c). 

FIG. 3: (Color online) (a) Configurations of the three-electron quantum dot with one spin- up 
electron occupying dot 1. (b) Three-electron energy levels calculated using the RSP-CI approach 
to the device shown in Fig. [T] (left-hand part) and using the appropriately fitted Hubbard model 
(right-hand part), (c) Similar spectra obtained with the LCAO-CI method 
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FIG. 4: (Color online) (a) Energy levels calculated by RSP-CI technique, measured from the 
ground state, as a function of the voltage applied to the control gate VG13. Black lines show 
energies of total-spin- 1/2 states, the red line shows the energy of the spin-3/2 state. Inset shows 
the three lowest energies as a function of the gate voltage, (b) Energies of three electrons localized 
on three Hubbard sites as a function of the tunneling amplitude t\3 measured from the ground 
state. Hubbard model parameters were extracted from exact diagonalization results. 

FIG. 5: (Color online) (a) Charging diagram of the triple quantum dot on resonance as a function 
of energy level E of each dot from N e = to N e = 6 electrons, (b) Stability diagram (E\ — [i = 
E2 — //; E3 — n) of the triple dot molecules with dots 1 and 2 tuned by a common gate. (N\, N2, N3) 
denotes average electron occupation and (— Ni, —N2, —N3) denotes average hole occupation of each 
dot. ji denotes the chemical potential of the leads. The light gray, yellow, and brown colors mark 
the stability regions of phases with total spin 0, 1/2, and 1, respectively. 

FIG. 6: (a) (Color online) Addition spectrum of a lateral triple quantum dot molecule measured 
by Gaudreau et al& Inset shows the layout of the gates defining the triple dot. (b) Elements of the 
triple-dot addition spectrum involving dots 1 and 2 only, drawn as a function of onsite energies of 
the two dots. The Hubbard parameters can be extracted directly from this diagram (see text for 
details). 



FIG. 7: (a) Addition spectrum of the lateral triple-dot device calculated within the Hubbard 
model after fitting to the experimental spectra shown in Fig. E^a). (b) (Color online) The same 
spectrum shown in the form of a charging diagram as a function of the single-dot onsite energies. 
Light gray, yellow, and brown color marks the stability region of molecules with total spin 0, 1/2, 
and 1, respectively. 
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